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We study the directional-ordering transition in the two-dimensional classical and quantum compass models 
on the square lattice by means of Monte Carlo simulations. An improved algorithm is presented which builds on 
the Wolff cluster algorithm in one-dimensional subspaces of the configuration space. This improvement allows 
us to study classical systems up to L = 512. Based on the new algorithm we give evidence for the presence of 
strongly anomalous scaling for periodic boundary conditions which is much worse than anticipated before. We 
propose and study alternative boundary conditions for the compass model which do not make use of extended 
configuration spaces and show that they completely remove the problem with finite-size scaling. In the last part, 
we apply these boundary conditions to the quantum problem and present a considerably improved estimate for 
the critical temperature which should be of interest for future studies on the compass model. Our investigation 
identifies a strong one-dimensional magnetic ordering tendency with a large correlation length as the cause of 
the unusual scaling and moreover allows for a precise quantification of the anomalous length scale involved. 

PACS numbers: 02.70.Ss, 05.70.Fh, 75.10.Jm 



I. INTRODUCTION 

The quantum compass model ^ has recently seen a re- 
naissance in condensed-matter physics, which was to a large 
part triggered by the observation that it may protect q-bits in 
a quantum computing setting |2, 3]. This observation may be 
of actual practical relevance as the quantum compass model 
can be realized by a special connection of Josephson junc- 
tion arrays, a concept with which first experimental successes 
could be reported A concrete realization in terms of real 
materials has also been proposed recently I5i|. Apart from 
the current interest from the quantum information perspec- 
tive, the quantum compass model is relevant as an effective 
description for orbital ordering, and was originally proposed 
in this setting |1|]. Due to the diverse interest in the model, 
recent contributions in the literature have studied many dif- 
ferent aspects, ranging primarily from detailed investigations 
of the ground-state properties 13, to a study of the possi- 
ble low-temperature phases and ^hase transitions in both the 
classical and quantum cases |7-9]. Complementary to that, 
recent studies considered modified variants of the compass 
model in one-dimensional chains lflol - fl2ll or in a magnetic 
field ifisll . In Ref. 13], two of us have proposed and stud- 
ied a two-dimensional (2D) geometric variant of the compass 
model. The model is also known to have relevance for other 
settings such as p + su perc onductors lfl5l[l6ll . the concept 
of dimensional reduction iflTll . and was recently shown to be 
isospectral 1 18] to Kitaev's toric code lfl9ll in a field. 

The 2D compass model (CJVl) is defined on a square lattice 
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of N ~ L X L sites as a (pseudo) spin model by the Hamilto- 
nian 

i i 

where Sf and are components of a two-component spin 
Si at site i. The spin can represent both classical and quantum 
degrees of freedom. In the latter case and are repre- 
sented by the usual Pauli matrices, i.e., S = (l/2)((Ta;, ay). 
The classical case is analogous to an ordinary classical XY 
spin S = {S'^.Sy) e The interesting feature of Eq. O is 
its anisotropy in spin and lattice space. 

For Jy ^ ,Jx, the ground state of Eq. ([T]i can be described 
by (weakly) coupled Ising spin chains oriented in the x or y- 
direction depending on | Ja;| > | Jy | or |Jy| > \Jx\, respec- 
tively. The quantum phase transition between these differ- 
entl y o riented ground states was shown to be of first-order 
1120112 ill . One interesting feature of that work is that Ref. ll2lll 
gives one of the first nontrivial applications of the recently in- 
troduced infinite pair-entangled tensor product states (iPEPS) 
ll22ll which aim at providing a new numerical approach to 2D 
interacting quantum systems. Following the same line of re- 
search, a quantum phase transition in a generalized CM has 
also been investigated recently ||23ll using the related multi- 
scale entangled renormalization ansatz (MERA) ll24ll . 

Here, we will focus on the symmetric case Jx ^ Jy ^ 
— 1 which allows - due to a discrete x ^ y symmetry in 
spin and lattice space ^ - for a thermal phase transition to 
a directionally-ordered low-temperature phase without long- 
range local magnetic order 10. Si- In Ref- 101, two of us have 
studied this transition extensively for both the classical and 
quantum CMs. One of the main results of this contribution is 
the confirmation that the CM suffers from extraordinary finite- 
size corrections when studied in a simple canonical ensemble 
on the torus, contradicting the naive assumption that generic 
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periodic boundary conditions are optimal. A solution to al- 
leviate this problem had already been suggested in Ref. 
where the authors proposed the use of so-called fluctuating 
boundary conditions (FBC) (sometimes also referred to as 
"annealed boundary conditions")- These formally place the 
CM in a larger configuration space where the partition func- 
tion is given by 

Zfbc — 

{Jb=±l} i 

instead of just 

ZPBC= flldS,cM-f3^CM). (3) 

i 

for the standard canonical ensemble with periodic boundary 
conditions (PBC). Here, { Jf,} denotes the set of boundary 
bonds on the periodic lattice which are allowed to fluctuate 
between —1 and +1 individually ll44ll . One assumes that the 
Jb degrees of freedom become unimportant in the thermody- 
namic limit. Indeed, it was shown in Refs. ^ and that 
FBC lead to very good finite-size scaling (FSS) properties in 
the classical case from which we have good evidence that the 
directional-ordering (DO) transition in the CM is in the 2D 
Ising universality class. 

Our good experience with these FBC is unfortunately of 
no use for the quantum CM because of the minus-sign prob- 
lem. Furthermore, simulations of the quantum CM are quite 
demanding and one may currently not reach large lattice sizes 
(say L > 64) with reasonable effort. In result, our current 
estimate of the critical ordering temperature is not very pre- 
cise as it rests on the use of non-optimal boundary conditions 
on moderate lattice sizes i^. Yet, given the large interest in 
the model we find it valuable to try to improve the accuracy 
of the critical temperature. A better knowledge of such quan- 
tities is necessary in order to tackle more advanced features 
such as the influence of disorder, etc. Jst]. Apart from the mo- 
tivation to improve the available critical data, there are further 
unsatisfactory points or open problems regarding the previous 
results 10, Ht] . These especially concern the ad hoc use of FBC 
to get precise results at the price of introducing extra degrees 
of freedom to the model. Why do these boundary conditions 
work so well and why do we observe a complete failure of the 
critical Binder parameter on periodic lattices? In this work 
we (re)address those questions and present improved results 
on critical properties of the classical and quantum CM that we 
obtained with a combination of algorithmic advances and by 
employing so-called screw-periodic boundary conditions. 

The outline of the rest of the paper is as follows. In Sec. |II] 
we start with a revision of our Monte Carlo (MC) approach 
and present an improved MC algorithm building on the Wolff 
cluster method. The latter will make possible a much more 
detailed comparison of FSS properties on periodic vs. fluc- 
tuating boundary conditions in Sec. |III1 going considerably 
beyond Ref. We will show that periodic boundary condi- 
tions behave even worse than previously anticipated. A solu- 
tion to this problem is thereafter suggested in form of screw- 
boundary conditions. These will allow to recover very good 



scaling properties without making use of an extended config- 
uration space (in form of fluctuating boundary conditions). 
Moreover they can be readily employed in quantum Monte 
Carlo (QMC) simulations which is the topic of Sec. lIVI where 
improved critical data for the quantum CM are presented. We 
end with a summary and conclusions in Sec. [V] 

II. OBSERVABLES AND MC APPROACH 

In this section, we present the standard approach to simulate 
the classical CM and describe in detail how we can improve 
the algorithm by making use of ideas from well-known cluster 
MC updates. A short discussion of the QMC approach for the 
quantum version of Eq. ([T]) is postponed to Sec. lIVI 

A. Revision of classical MC approach and relevant observables 

In the classical case, we have investigated the ensembles 
specified by Eqs. (|2]i and (O using the Metropolis algorithm 
combined with the parallel tempering (PT) scheme 11251 - 12711 
parallelizing simulations at different temperatures. Technical 
details of this approach are described in Ref. J^. During a MC 
simulation, we measure an order-parameter known to describe 
directional-ordering 

D^j^\E,-Ey\, (4) 

with = Sf+e. and Ey = ^f+e„- We concen- 

trate here on its susceptibility 

X = N {(D^) - (Df) , (5) 

which diverges at the phase transition temperature Tc- On 
finite systems the divergence in x is smoothened into a fi- 
nite maximum Xmax(i) at some pseudocritical temperature 
T„iiix{L). Finite-size scaling predicts the following two fun- 
damental scaling relations (see, e.g., Refs. ll28l - t30ll ) 

Xmax(i) - L-y/-" , (6) 
T,nax(i) = Te + aL-^/'' , (7) 

which are the primary means used in this paper to obtain the 
critical exponents i', 7 and the critical temperature Tc, and to 
discuss anomalous scaling. From MC simulations at discrete 
temperatures in the vicinity of the phase transition, we obtain 
Xmax(-^) and T„iax(i) by making use of standard reweighting 
techniques llsill and optimization algorithms. Error estimates 
for these quantities are obtained by "jackknifing" this proce- 
dure ||30il32ll. 

B. A Wolff-like cluster update 

Up to now the application of the PT technique has proven to 
be quite efficient, enabling a study of the classical (and quan- 
tum) CM on moderately large lattice sizes i^. However, for 
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linear system sizes of about L > 100, we observe that the 
method runs into problems as the equilibration times in the 
MC simulation become visibly very long. In order to go effi- 
ciently beyond such lattice sizes, a further improved method 
is therefore called for. 

Indeed, as we shall propose here, a rather straightfor- 
ward improvement is possible with a special one-dimensional 
Wolff-cluster update [33]. To see this reconsider the ordinary 
Wolff-construction for 0{N) spin models with Hamiltonian 
'Ho{N) = JJ2{ij) ^i^j- Following Ref. ill], the operation 



R'^iS,) = S, - 2(S, • r)r 



(8) 



denotes a reflection of the [0{N)] spin along a hyperplane 
defined by the vector r. Given a random r, Wolff clusters are 
constructed using the bond activation probability 

P,, = 1 - cxp (min[0, - J/3[S,S, - S,i?'^(S,)]) (9) 

for bonds (ij). The spins in each cluster are then flipped by 
applying Si — > R'^{Si) which implements the (non-local) MC 
move. The principle of detailed balance is satisfied by requir- 
ing the invariance 



iJ(i?"-(S,),i?'-(S,)) = H{S,,Sj) = H,. 



(10) 



of the bond energy Hij (H — ^ij) under reflection of 

the spins for each bond (ij) of the lattice. While this is true 
for 'Ho{N)^ it is clearly not true for the CM in general. How- 
ever, we know that for the CM, the following special reflection 
operations i?*'^ and with 



iSf, -Sf) 



(11) 

(12) 



are symmetries on a subset of all bonds 12,1^, namely 
fliat 
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FIG. 1 : (Color online) (a) Comparison of the autocorrelation time r 
at the critical point using periodic boundary conditions. Comparison 
of t(L) for the pure Metropolis and the combined Metropolis + clus- 
ter update. An overall improvement for r as well as better scaling is 
evident for the cluster variant, (b) Comparison of r for L = 36 as ob- 
tained from the pure parallel tempering approach with the improved 
parallel tempering variant. 



H{R-y{S,),R-y{S,+,J) = H{S,,S,+,J, (13) 
l/(i?^==(S,),i?"^(S,+eJ) = ff(S„S,+eJ. (14) 

Thus, we may employ R'^" and R'^" to construct one- 
dimensional clusters of spins along the x- or y-direction. Em- 
ploying the form of the CM Hamiltonian ([TJ and the general 
relation (|9]l we obtain the following bond-activation probabil- 
ities 

P,,+e, = 1 - cxp(-2J,/35f5,%J, (15) 
Pu+e, = 1 - cM-^JyPSfSf^J (16) 

for cluster growth along the x- and y-directions, respec- 
tively. Note that the cluster construction is really strictly one- 
dimensional, i.e., when we build x-clusters we do not attempt 
to add y-bonds to the cluster which would break condition 
(fT3T i. Cluster construction starts as usual by picking a random 
start site from which cluster growth proceeds. 

An obvious difference to the original Wollf algorithm is the 
discrete set of possible spin reflections. Thus the cluster up- 
date alone does not satisfy ergodicity. This is not a problem as 



long as ordinary Metropolis (as well as PT updates) are per- 
formed in addition. In each MC sweep we perform on average 
L cluster updates in both x and y directions as well as N local 
Metropolis updates. We have verified by detailed comparison 
to existing data that the new algorithm works correctly. Let 
us proceed directly to an evaluation of the new update. In or- 
der to examine its performance we ran several tests on lattice 
sizes L = 16, 24, 32, 48, 64, (96 in case of the cluster update) 
in the ensemble ZpBc at the pseudocritical temperatures 
Tmax{L) (known from our previous study). In the first test 
we switched off the PT update and compared the autocorrela- 
tion time r of the energy time-series which should scale at the 
critical point like t ^ L^. Figure[TJa) compares the scaling of 
T with and without the above cluster updates. Clearly, we find 
that the cluster algorithm behaves much better Apart from 
the expected absolute reduction of r we observe a decrease 
of z from z w 3.5 to z w 2.5 which is apparent from the 
different slopes in the log-log plot. Next, we also compared 
the autocorrelation time t in simulations employing the PT al- 
gorithm. Without performing a detailed scaling analysis, it is 
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evident from Fig.lHb) that the cluster update further improves 
the PT algorithm. 

The methodological improvement presented here allows to 
study much larger system sizes than before. In the course of 
this study, we have performed simulations up to L = 512. 



III. CLASSICAL COMPASS MODEL: RESULTS 

In this section we are going to employ the algorithmic ad- 
vances to restudy critical properties of the classical CM. Spe- 
cial focus is given to a more detailed comparison of ensembles 
Zfbc and ZpBC- Based on this comparison we will thereafter 
propose the use of alternative boundary conditions and study 
their effects on FSS properties. 



A. Revisiting periodic and fluctuating boundary conditions 

Previous investigations of the DO transition have clearly 
shown that the use of ensemble ZpBC is favorable over ZpBc 
in terms of FSS properties 10,01, where the most severe "fail- 
ure" of ZpBC establishes itself in an unconventional behavior 
of the Binder parameter. Despite these observed defects it was 
argued 121] that one may still employ PBC to extract the criti- 
cal properties given the system sizes L are large enough. This 
argument was supported from extrapolations of pseudocritical 
temperatures T,„ax(i) which gave consistent results for both 
Zfbc and Zpbc of about = 0.1464(2). 

With the newly available cluster procedure, we will inves- 
tigate this issue further to make more quantitative statements 
about how ensembles ZpBC and Zpbc converge towards an- 
other asymptotically. We have thus simulated the CM for sys- 
tem sizes between L = 12 and L = 512, pushing L a factor 
of 5 — 10 times larger than before. In comparison to Ref. i^, 
we have added system sizes L ~ 96, 128, 164, 256, 512. 

The observables described in Sec. Ill Al were estimated us- 
ing about 10^ samples. We have taken measurements only 
every m MC sweep such that the final autocorrelation time 
was small, r < 10 (to was in the range of 4 — 100). For 
a presentation of the typical temperature dependence of the 
order parameter and susceptibility we refer the reader to the 
previous work of Ref. i^. Here we just present the pertinent 
data obtained for the pseudocritical temperature Tmax(i) and 
Xmax(i)- Figure |2] summarizes the FSS analysis for the two 
different ensembles considered. The partly surprising results 
of this comparison can be summarized as follows. 

First, we observe that the FSS behavior for Zfbc is fully 
consistent with earlier results, i.e., data obtained on larger lat- 
tices agree with the extrapolations from smaller lattice sizes. 
This further confirms the claim of 2D Ising universality be- 
yond any reasonable doubt. Indeed, fits to Eq. (|7]i yield our 
estimate of the critical temperature and critical exponent as 



Tc = 0.14621(2), 



1.00(1) 



(17) 



(with x^/d.o.f = 1.4) which represents an improvement of 
roughly one order of magnitude over the previous estimate. 
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FIG. 2: (Color online) Improved finite-size data for the classical 
compass model with periodic (PBC) and fluctuating boundary con- 
ditions (FBC). (a) Extrapolation of the pseudocritical temperatures 
Tmax(i'). Data for FBC follow a perfect straight line. The periodic 
case shows a clear bend for system sizes L > 96 not anticipated 
before based on the straight line extrapolation in Ref. (b) The 
susceptibility maxima divided by A'^ behave non-monotonously and 
indicate a resonance phenomenon at about L ~ 120 (indicated by 
the arrow), (c) FSS of Xmax in a log-log plot. For FBC a power law 
is evident with ^/u = 1.75(1). For PBC a power-law extrapolation 
is not justified and different scaling regimes are apparent. 
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Together with the critical exponent 7 = 1.75(1) obtained 
from analyzing the scaling of Xmax^ this is in perfect agree- 
ment with the exactly known critical exponents for the 2D 
Ising model. 

Second — and this is the surprising result — the scaling 
for ZpBc reveals a more complicated or stronger anomalous 
scaling than previously thought. This is especially apparent 
in Fig. 121 a) where pseudocritical temperatures for L > 96 
clearly deviate systematically from the previous extrapolation 
[upper (blue) straight line] of Ref. |@] based on the assumption 
of 2D Ising scaling for L > 30. Note that the upper (blue) line 
was also justified because it matched exactly with the result 
from ZpQc- However, instead of a clear power law scaling 
in 1/L we observe a "double bend" in the FSS curve which 
seems to collapse onto the curve from Zfbc for very large 
system sizes L > 256 (inset of Fig. |2 a)). 

The same anomalous scaling behavior shows up in the 
quantity Xma,x{L)/N of Fig. |2b) which shows a strong non- 
monotonic behavior at a length scale of about L « 120. Thus, 
any attempt to extract the critical exponent 7 from a log-log 
plot as in Fig. |2c) is doomed to fail on length scales below 
L « 256. The observed non-monotonic behavior also shows 
up in the Binder parameter but we postpone a discussion on 
that to the next subsection. 

As a matter of fact, it is thus totally unreliable to obtain crit- 
ical exponents and the critical ordering temperature Tc from 
simple extrapolations in the ensemble Zpec on periodic lat- 
tices (at least for L < 256). The previous seemingly cor- 
rect extrapolation was a matter of coincidence. Turning this 
observation around, one might even be tempted to argue for 
non-Ising behavior in the CM if one did not have access to the 
largest lattice sizes studied here. This situation is most unsat- 
isfying and calls for a deeper investigation and a workaround. 
We will attempt precisely this in the next subsection. 



B. Screw-periodic boundary conditions 

The main message of the discussion so far is that PBC show 
a more complex scaling behavior than previously thought with 
the appearance of a clear resonance effect and at least two dif- 
ferent scaling regimes. This disqualifies the use of PBC to 
extract critical properties. The FBC ensemble on the other 
hand also has — despite its intriguing performance — a cou- 
ple of drawbacks. The most important of all is that we may 
not easily use it in QMC because fluctuating couplings induce 
a minus-sign problem. Second, one may wonder whether it is 
safe to use them in the first place as a trustable FSS theory is 
not available and one is actually simulating a different model. 

Here, we would like to ask whether it is nevertheless pos- 
sible to deal with the described problem using only slightly 
modified boundary conditions without going to a higher di- 
mensional configuration space. 

It is quite obvious that the torus geometry hides or shields 
the true physics going on. One possibility to unveil the true 
properties of the CM in the thermodynamic limit is to intro- 
duce systematic deformations to the torus. In this way we can 
at least see how the problem is alleviated (or made worse). 
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FIG. 3: (Color online) Illustration of screw-periodic boundary con- 
ditions along the y axis as defined in Eq. JI8t . Two examples with 
S — 1 and 5* = 2 are shown. In our simulations the same procedure 
is applied to the a::-direction. 



Among all such deformations one may consider a Mobius 
strip or so-called screw-periodic boundary conditions. Such 
deformations are very easy to implement on the computer 
and cost no extra updates. We decide to study screw-periodic 
boundary conditions (SBC) which are defined by 



{x,y + l) 



{x,y+l) \f y < L 

{[x + S] modL,0) if y = L 

(x + 1, y) if X < L 

(0, [y + S] mod L) if x^L 



(18) 



where, e.g., {x,y + 1) denotes the nearest-neighbor of site 
i = {x,y) in y-direction. The parameter is a parameter 
that determines how much we deform the clean torus case. 
Figure [3] illustrates this concept for two cases S ~ 1 and 
S = 2. The cases 5 = and S* = L are obviously iden- 
tical to the usual PBC. For a given lattice size, S may take 
only certain values in order to satisfy the overall periodicity 
constraint. The possible S'-values are given by the set of all 
(distinct) divisors of L. SBC are discussed in various forms in 
the literature, sometimes called helical boundary conditions or 
shift boundary conditions. Mostly they have been employed 
because theyhave some advantages regarding implementa- 
tion issues II29I [34I1 or to complement FSS analysis Issll as 
they approach the thermodynamic limit with (slightly) differ- 
ent pseudocritical temperatures. A further useful application 
is the controlled formation of tilted domain walls in the Ising 
model 136]. Note that each site still has exactly four neighbors 
which distinguishes SBC from open boundary conditions. 

SBC allow one to put the lattice points into representation 
classes which we will call loops. A loop is the set of all points 
i that the screw/helix passes until it closes itself. The length 
of a loop is called Li and is the number of points it contains. 
Each point i is obviously member in exactly two loops, one 
for the X- and one for the y-direction. Given a lattice size 
L and compatible screw parameter S, each loop has length 
Li = L"^ / S and for S = 1 all points belong to the same 
loop. The notation introduced here becomes relevant when 
discussing the symmetry properties of the CM and its ground- 
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FIG. 4: (Color online) Dependency of the susceptibility xiT) for 
L = 36 on the choice of the screw parameter S. The case S* = 36 
is equivalent to periodic boundary conditions. With decreasing S 
(or increasing the boundary loop length) a clear shift in the peak 
is observed with an apparent resonance at S ~ 12. For 5=1 the 
susceptibility is (nearly) identical to the susceptibility obtained in the 
fluctuating bond ensemble Zybc (continuous line without symbols). 
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FIG. 5: (Color online) Finite-size behavior of the Binder parame- 
ter B at the critical point T — Tc for periodic (S — 0) and screw 
boundary conditions with S = 1. While there is a strong anomaly 
for S — 0, the case 5=1 almost completely removes the defect and 
is consistent with a scale independent value of B at the critical point. 
Moreover, it is consistent with the usual 2D Ising value (jTllSal in- 
dicated by the horizontal line. Note that the symbol size for 5 = 1 
is bigger than the error bar. 



state degeneracies. 

A simple check confirms that the usual one-dimensional 
spin flip operators P/ ^ J]^ crf^ i) ™d Qm = Uy f^f,„ y) 
H il (which are related to operations R'^'^ and R'^y in 
Sec. Ill Bl i are no longer symmetries of the (quantum) Hamil- 
tonian ^ if S ^ 0. However they can be generalized to the 
SBC case with the following operators 

Pi =11-1 (19) 

Qm =l[^f, (20) 

where I and m now refer to a loop along the x or y-direction. 
As we can control the number of independent loops via the 
parameter S, we can control the number of such symmetry 
operators and thus the degeneracy of the ground-state. Indeed 
it is possible to change the ground-state degeneracy from ex- 
ponential growth 2^+^ {S = 0) to a constant 2{S = 1), an ob- 
servation which may have interesting physical consequences. 
The thermal DO transition studied here should not be affected 
by this as the relevant global Z2 symmetry is not changed. 

Let us proceed to study the effect of SBC in actual MC sim- 
ulations. To this end, we choose a system size L = 36 which 
allows us to study quite a large number of screw parameters 
5 = 0, 1, 2, 3, 4, 6, 9, 12, 18. In each case we measured the or- 
der parameter D and its susceptibility x for a couple of tem- 
peratures close to the phase transition. Figure |4] depicts the 
drastic effect of SBC on the susceptibility x- Starting from 
the periodic case S' = (or S* = 36) we observe that x moves 
massively towards the curve from Zfbc for decreasing 5 > 
or increasing Li. The case S* = 1 gives an almost identical 
result to that obtained with fluctuating couplings in the Zfbc 
ensemble. Second it is apparent that there is a resonance at 



some length scale determined by S* « 12 at which the fluctu- 
ations in the system are strongest. 

The above picture thus suggests that S* = 1 resolves the 
FSS problems observed in the CM for the susceptibility al- 
most completely. Furthermore it gives a hint at the order of 
the disturbing length scale (L; ?» 36^/12 « 110) which is 
present and which prohibits the extraction of correct critical 
data. Any solution that restores good FSS properties should 
also repair the behavior of the Binder parameter 

1 (D^) 

^^^-3(W' ^''^ 

whose normally used power is due to a scale invariance at 
the critical point with only leading order corrections. Thus, 
if SBC really solve the problem they should also remove 
the very unconventional behavior in the Binder parameter 
which was observed with PBC in Ref. [9]. Figure |5] shows 
a comparison of the finite-size behavior of B for the cases 
S ~ (periodic) and S = 1 performed close to the criti- 
cal point given in Eq. ( fTTj i. The periodic case shows the ex- 
pected non-monotonous behavior (with a possible restoration 
for L > 256). The 5* = 1 screw restores the expected scal- 
ing behavior - up to a small bump for L < 12 - completely, 
i.e., it is almost a constant for various system sizes and agrees 
rather well with the known value of B « 0.61 IstIIssIi for the 
2D Ising model (constant line in Fig.|5]). Note, however, that 
the agreement is not expected to be perfect as boundary con- 
ditions can have (a small) influence on the (only weakly uni- 
versal) critical value of B ll39ll . For an analysis of the Binder 
parameter for FBC we refer the reader to Ref. 

This brings us into the position to claim that SBC are a very 
efficient tool to study critical properties of the CM. Before we 
apply these to the quantum CM let us try to shed some light 
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onto the origin of anomalous scaling (with PBC). 



Origin of anomalous scaling: one-dimensional spin order 



It is evident from the MC analysis in Sees. 1111 Al and HITB] 
that there is a second length scale in the CM which influences 
fluctuations and which can be overcome by SBC. Let us now 
turn to a discussion of possible reasons for this as a more fun- 
damental understanding of this phenomenon is clearly desir- 
able. 

We know that the low-T (directionally-ordered phase) of 
the CM is essentially one-dimensional where the spins along 
each row or column are essentially decoupled. These spins 
thus form a ID spin chain. Using this picture, a plausible 
explanation for the failure of FSS was actually suggested in 
Ref. 101 where it is argued that the magnetic spin-spin corre- 
lation length along each chain exceeds the linear system 
size L at low temperatures. If this were the case, all spins 
would align themselves along each chain although a direc- 
tionally ordered state can survive even with domain walls in 
spin space. Such total magnetic ordering tendency could in- 
fluence the fluctuations of the true order parameter, making it 
more robust against thermal fluctuations and spoiling its FSS 
properties. 

To test this hypothesis let us write down an order parameter 
for such one-dimensional magnetic ordering tendency Mm 

as 



Ml 



D 



x(y) — loops,/ 



1 



(22) 
(23) 



Here Nl denotes the number of boundary loops (as intro- 
duced above) and Li the length (number of sites) of loop /, 
i.e., we already take care for the general screw-periodic case. 
The quantity 1\Iid probes whether all spins along each chain 
(or loop) like to align themselves. 

To test whether such possible (long-range) ordering of the 
spins exists on top of directional order, we perform a couple of 
MC runs at the critical temperature obtained in Eq. ( fTTI ) for 
lattice sizes L = 8, 12, 20, 36, 48. In each case we simulate 
all possible screw parameters S. In Fig.|6]we plot the expecta- 
tion values of Mio vs the screw loop length L; = I? j S. Re- 
markably, the data from different system sizes collapse onto 
the same curve for Li < 100 where a finite expectation value 
for M\Y) is evident. This magnetic ordering tendency does 
not persist in the thermodynamic limit as for L; > 100 it 
suddenly approaches 0. We conclude that there is a strong 
tendency for the spins to align themselves which is enforced 
by PBC. Application of SBC can overcome this problem be- 
cause it exceeds the typical length scale along each loop. The 
same is true for FBC by artificially introducing kinks in the 
spin configurations, which is the basic reason why they do not 
show anomalies such as those in Fig. 2(b). 

The length scale at which the sudden decrease in Aim 
appears coincides precisely with the non-monoticities ob- 
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FIG. 6: (Color online) Expectation values for one-dimensional mag- 
netization for several lattice sizes and choices of the screw displace- 
ment S at the critical temperature Tc . The x-axis is the loop length 
Li = /S. Data from different lattice sizes collapse onto one curve 
(for Li < 100). A clear crossover from a region with finite mag- 
netization to a disordered spin state is observed on a length scale 
Li ~ 100 (indicated by the arrow). This length scale corresponds to 
those where anomalies are seen in the FSS analysis. 



served in the scaling of Tmax and Xmax of Fig. |2] The res- 
onance effect in Fig.|4]can be explained because at Li = 
we have strong fluctuations in Mid in addition to the normal 
fluctuations in the directional-order parameter D. 

These results essentially confirm the picture of Ref. |01 and 
quantify precisely the length scale involved. The quite large 
magnetic correlation length can be understood by recalling the 
exponential divergence of the magnetic correlation length at 
low temperatures in the ID Ising model. 



IV. RESULTS FOR THE QUANTUM CASE 

We have now developed everything to proceed to the main 
objective of this paper which is to improve the estimate of 
the critical ordering temperature Tc for the DO transition in 
the presence of quantum fluctuations. Due to the results of 
Sec. 1111 Al it is probable that the previous result Tc = 0.055(1) 
in Ref. tg] is slightly off the true critical temperature due to the 
presence of the second length scale. 

We expect that SBC rectify and improve this value. There- 
fore, new QMC simulations in the stochastic series expansion 
(SSE) framework using directed-loops ll40l liHl and FT up- 
dates were performed implementing S* = 1 SBC. Otherwise 
our approach rests on that presented in Ref. |9|] where concrete 
implementation issues are discussed. A couple of simulations 
for lattice sizes L = 10, 12, 16, 20, 24, 28, 32, 42 were per- 
formed and approximately 100 000 statistically independent 
samples of the order parameter D were taken in each case. 
The pseudocritical temperatures Tmax(-^) were obtained from 
the peak in the variance of D utilizing the quantum general- 
ization of the multihistogram reweighting idea ll42ll . Figure 
|7] shows the pseudocritical temperatures obtained and com- 
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FIG. 7: (Color online) FSS plot of the pseudocritical temperatures for 
the quantum CM from the susceptibility comparing periodic (S = 0) 
and screw-periodic boundary conditions (S = 1). The latter clearly 
lead to a considerable improvement. The line is a fit using a power 
law correction in l/L. The dashed extrapolation of data from peri- 
odic boundary conditions underestimates the critical temperature 
as expected from the discussion in Sec. lIII Al 
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no. 8: (Color online) Analysis of the one-dimensional magnetiza- 
tion M^jj for the quantum CM. The result from the classical case is 
indicated by the line taken from Fig. [6] (and divided by a factor 2). 
The magnetic length scale is clearly much larger in the classical case. 



pares them to the old data utilizing PBC. As expected, an 
evidently improved FSS behavior is observed for the screw- 
periodic case. This is apparent from the absolute move of 
Tmax towards the true critical temperature for small L and the 
much better power-law scaling in terms of l/L. Indeed, the 
SBC data are fully consistent with v — 1 and a straight line fit 
to 

T^UL)=T, + aL-^ (24) 

yields our new estimate for the critical temperature as 

Tc = 0.0585(3) (25) 

with x^/d.o.f = 1.5 using all lattice sizes studied. Even 
using only the of smallest systems L = 12 to L = 20 the 
extrapolation yields a consistent value of Tc = 0.058(1), a 
property which is of most practical relevance for studies aim- 
ing at numerically verifying more qualitative effects (see, e.g., 
Ref. i^). Leaving as a free fit parameter as in Eq. (|7]i we 
obtain Tc = 0.0586(8) and v = 0.97(15) which is consis- 
tent with 2D Ising behavior. Hence, although we have per- 
formed much less simulations and on smaller system sizes, 
we have obtained a much better and improved result just by 
an adequate choice of the boundary conditions. The present 
result does not agree within error bars with the previous es- 
timate Tc ~ 0.055(1) because of the anomalous behavior 
which was not accounted for in the ordinary Ising extrapola- 
tion (dashed line in Fig. 13 with a correction on periodic 
lattices. However, it appears that the effect of the magnetic 
length scale is not as severe as in the classical case. This could 
be expected due to the presence of quantum fluctuations. On 
the other hand the temperature regime is lower which could 
in principle even stabilize the unwanted order. In order to 
get an approximate estimate for the length scale involved, we 



have finally analyzed the one-dimensional magnetization also 
for the quantum case, where we restrict ourselves to measure 
^Iijj (which corresponds to the quantization direction) along 
i/-loops. Figure [8] shows the SSE estimates for M^jj for var- 
ious system sizes and screw parameters at T = 0.07 (chosen 
for convenience because it is close to but still in the region 
where PBC show unusual behavior). It verifies that quantum 
fluctuations reduce the overall value of M^jj and that they 
lead to a clear diminution of pseudo-magnetic order at a scale 
corresponding to roughly L « 50 which is apparently smaller 
than in the classical case (line in Fig. [S]). Moreover, we also 
arrive at this conclusion by studying the behavior of the sus- 
ceptibility (similar to that in Fig.|4]i and find that the resonance 
is shifted to a smaller length scale in accordance with the find- 
ings just described. However, even such a moderate scale can 
still be a formidable challenge to overcome for QMC without 
SBC. 

In summary, the estimate of Tc for the quantum CM has 
seen several steps of adjustments on a relatively short time 
scale as summarized in Table U The result of this section pro- 
vides the new benchmark estimate which should be useful for 
future studies. 



TABLE 1: Previous and current estimates of the critical temperature 
of the DO transition exemplifying the (previous) difficulty of its ex- 



traction. 


Tc 


System sizes 


Boundary cond. 


Method 


Ref. 


0.075(2) 


10 


-20 


Periodic 


Trotter QMC 


La 


0.055(1) 


10 


-96 


Periodic 


SSE + PT 




0.058(1) 


12 


-20 


Screw 


SSE -1- PT 


This work 


0.0585(3) 


12 


- 42 


Screw 


SSE -1- PT 


This work 
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V. SUMMARY AND CONCLUSIONS 

Summarizing, we have revisited the directional-ordering 
transition in the classical and quantum compass models em- 
ploying two types of methodological advances. 

In the classical case we were able to formulate a spe- 
cial one-dimensional cluster update which in combination 
with Metropolis and PT methods allowed to investigate much 
larger system sizes than before. The following detailed com- 
parison between the classical CM with periodic boundary con- 
ditions and a fluctuating bond ensemble showed that periodic 
boundary conditions scale much worse than known so far In- 
stead of the usual power law, anomalous scaling becomes evi- 
dent with a resonance and non-monotonic behavior in the sus- 
ceptibility and the Binder parameter at length scale of about 
L ?» 100 — 200. In any typical MC simulation one would 
therefore not be able to predict critical properties correctly 
when the simulation is done with periodic boundary condi- 
tions. This resonance is argued to be due to a magnetic corre- 
lation length which prohibits the formation of domain walls 
at finite temperature on small clusters. To counteract this 
problem we have proposed to employ screw-periodic bound- 
ary conditions. We have shown that they are able to remove 
scaling anomalies in the classical case almost completely. 

This concept then proved to be a key step for simulations 
of the quantum compass model where we were able to obtain 



a more accurate estimate of the critical DO temperature based 
only on the change in boundary conditions. On the physical 
side we have seen that the CM represents a formidable chal- 
lenge despite its simplicity — even for well settled numerical 
approaches. The right choice of boundary conditions or topol- 
ogy is more essential for numerical studies of the CM than for 
many other models. 

Technically, it is clear that screw-periodic boundary condi- 
tions should be used in future studies of various other aspects 
in the quantum compass model. Moreover, we regard SBC as 
a well suited and general method which deserves more atten- 
tion even in studies of other systems. Via the screw parameter 
S one may be able to tune or minimize corrections to scaling. 
We are currently applying them to further studies of the quan- 
tum phase transition in 2D dimerized Heisenberg models (see, 
e.g., Ref. iH). 
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